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Recently, we have proposed a novel computation method of first-passage times of a ran- 
dom walker between a starting site and a target site of regular bounded lattices of arbitrary 
shape P]. Such first passage time properties are for example crucial to describe the kinet- 
ics of diffusion limited reactions in confined media . The obtained expressions involve 
pseudo-Green functions[4], which have been estimated according to different approximation 
schemes. In this note, we give the exact expression of these pseudo-Green functions in both 
cases of a rectangular domain with reflecting boundaries and a rectangular domain with 
periodic boundary conditions. This allows us to provide exact and explicit expressions of 
mean first-passage times with one or two target sites. 

Consider first a random walker starting at the source S, of position of a rectangular 
regular lattice with reflecting boundary conditions. The mean time it takes to reach the 
target T, of position r^, for the first time, is given by[l| 

{T) = N[H{rT\rT)-H{rT\rs)] (1) 
where H is the pseudo-Green function, which satisfies 

{k,i) 

a being the coordination number of the lattice, namely 4 for a 2-D lattice, or 6 for a 3-D 
lattice. In this expression, the sum runs over all the neighbours k of the site i, considering 
that a site which is at the boundary of the domain is its own neighbour. is the total 
number of sites of the lattice. In the presence of two absorbing targets, the mean time it 
takes to reach either of the two targets is: 

(T) — j\r (-^oi ~ His){Ho2 - H2s) - (-^12 - -^25) (-^12 - Hu) 

Hqi + Hq2 — 2i/i2 

while the eventual hitting probabilities Pi to reach the target i writes 



Hqi+Ho2-2Hi2 

H2s+Hoi-Hi,-Hi2 

Ho'L+Ho2-2Hx2 



(4) 



where H12 = H{rT^\rT2) and, for i = 1 or 2, His = H{rT^\rs), H^i = H^Tt^IttJ. 

The exact expression of the pseudo-Green function H involved in previous equations may 
be computed explicitely with the help of Fourier analysis. For a 2D domain with X sites in 
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FIG. 1: 2D Two-target simulations (red crosses) vs. theory (plain line). One target is fixed at 
(-5,0); The source is fixed at (5,0); The other target is at (x,3). The domain is a square of side 41, 
the middle is the point (0,0). The absorption time is normalized by the number of sites A^. 



the X direction and Y sites in the y direction, H writes 
H{r\r') 
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where x and y are the coordinates of r and x' and ?/' those of r'. Here, the coordinates of the 
left-bottom corner are taken equal to (1/2, 1/2), so that all the coordinates are half-integers. 

Figure n] shows the mean absorption time (T) given by eq. Q as a function of the target 
positions together with numerical simulations in the case of a 2D square domain. These 
simulations have been performed with a method based on the exact enumeration methodj^. 

Similar results can be obtained in the case of a rectangular domain with periodic boundary 
conditions. Here, the pseudo-Green function is 

T X—\ Y—1 2irmT(x—x') 2niTT(y—y') 

H(r\r') = — V V X Y 

^^^^ ' N ^ 1 _ 1 (cos^ + cos^ 
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where 5om is the Kronecker delta. 

Note that our results can be easily extended to the 3D case. 

We gratefully thank S. Redner for suggesting us the simulation method. 
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